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ABSTRACT 

We report ultraviolet spectra of the high-redshift (z cm « 2.9) quasar, HE 2347—4342, taken by the 
Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST). Spectra in the G130M 
(medium-resolution, 1135-1440 A) and G140L (low-resolution, 1030-2000 A) gratings exhibit patchy 
Gunn-Peterson absorption in the 303.78 A Lya line of He II between z = 2.39-2.87 (G140L) and 
z = 2.74-2.90 (G130M). With COS, we obtain better spectral resolution, higher-S/N, and better 
determined backgrounds than previous studies, with sensitivity to abundance fractions XHeii ~ 0-01 
in filaments of the cosmic web. The He II optical depths from COS are higher than those with the 
Far Ultraviolet Spectroscopic Explorer (FUSE) and range from th c ii < 0.02 to th c ii > 5, with a 
slow recovery in mean optical depth to (th ii) < 2 at z < 2.7. The He II/H I optical-depth ratio 
varies (17 « 10-100 for 2.4 < z < 2.73, and -q = 5-500 for 2.75 < z < 2.89) on scales Az < 0.01 
(10.8 Mpc in comoving radial distance at z = 2.8), with numerous flux-transmission windows between 
1135-1186 A. The He II absorption extends to 1186.26 A (z = 2.905), including associated absorbers 
with z a b s « zqso and minimal "proximity effect" of flux transmission at the He II edge. We propose 
a QSO systemic redshift zqso = 2.904 ± 0.002, some Az = 0.019 higher than that derived from 
O I A1302 emission. Three long troughs (4-10 A or 25-60 Mpc comoving distance) of strong He II 
absorption between z — 2.75 — 2.90 are uncharacteristic of the intergalactic medium if He II reionized 
at z r « 3. Contrary to recent indirect estimates (z r — 3.2 ± 0.2) from H I optical depths, the epoch 
of He II reionization may extend to z < 2.7. 

Subject headings: quasars: reionization — individual: HE 2347-4342 galaxies: active, intergalactic 
medium, quasars: absorption lines, ultraviolet: general 



1. INTRODUCTION 

The epoch of reionization in hydrogen has become a 
topic of considerable interest (Barkana & Loeb 2001; Fan, 
Carilli, & Keating 2006; Meiksin 2009; Furlanetto et al. 
2009) as a probe of the transition from neutral to ionized 
hydrogen in the intergalactic medium (IGM) . This tran- 
sition occurred somewhere between redshifts z = 6-12, 
marking the exit from the cosmic "dark ages" , beginning 
at the time when the first stars and galaxies formed at 
redshifts z > 30 (Tegmark et al. 1997; Ricotti, Gnedin, 
& Shull 2002; Trenti & Shull 2010). Helium underwent 
similar reionization from He II to He III (that is, from 
He+ to He+ 2 ) at z = 2.8 ± 0.2 (Reimers et al. 1997; 
Shull et al. 2004), most likely mediated by the harder 
(E > 54.4 eV) radiation from quasars and other active 
galactic nuclei (AGN). With a 4 ryd ionization potential, 
Hc + is harder to ionize than H°, and Hc +2 rccombincs 
5-6 times faster than H + (Osterbrock & Ferland 2006; 
Fardal, Giroux, & Shull 1998). For these reasons, and 
the fact that most hot stars lack strong 4 ryd continua, 
it is believed that AGN are the primary agents of He II 
reionization. Ionization models find that, owing to its 
resilience, He + is much more abundant than H , with 
predicted column-density ratios r\ = N(He II) /N(H I) « 



50-100 (Miralda-Escude et al. 1996; Fardal, Giroux, k 
Shull 1998). 

Determining when and how the universe was ionized 
has been an important question in cosmology for decades 
(Gunn & Peterson 1965; Sunyaev 1977). Although recent 
progress has led to quantitative constraints, we still do 
not know whether galaxies are the sole agents of hydro- 
gen reionization, and the epoch of reionization remains 
uncertain. We do know that hydrogen reionization in 
the IGM is complete by z « 6, based on the rapid evolu- 
tion between redshifts z — 6.2 and z ~ 5 of Lya Gunn- 
Peterson (GP) absorption from neutral hydrogen along 
lines of sight to QSOs (Fan et al. 2002; Gnedin 2004). 
The optical depth of the cosmic microwave background 
(Komatsu et al. 2010) sets a 3er constraint on instan- 
taneous reionization at z r > 6.5 and a la confidence 
interval of z r = 10.5 ± 1.2. A limit z r > 7 is set by the 
detection of Lya emitters at z — 6.5 (major evolution 
in the observed luminosity function would be expected 
if neutral gas were present around the emitters). These 
observations are in line with ACDM simulations (Trac & 
Cen 2007; Cen 2003; Gnedin 2004) and ionization-front 
(I-front) models (Venkatesan, Tumlinson & Shull 2003; 
Bolton & Haehnelt 2007; Shull & Venkatesan 2008). The 
current data leave open several hydrogen ionization see- 
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narios, some involving simple reionization at z r w 10 
and others with more complex ionization histories that 
depend on the star-formation rate at z = 7 — 20 and the 
transition from Pop III to Pop II stars. Numerical mod- 
els of the reionization process often follow the expansion 
and overlap of I-fronts, but the details are complicated 
by geometric uncertainties in the ionizing sources, higher 
recombination rates in denser gas, and radiative transfer 
of Lyman continuum (LyC) photons. In particular, the 
gas is highly structured in both density and temperature. 

The He II reionization epoch (z r ) and ionization pro- 
cesses are also poorly characterized. The thermal evolu- 
tion of the IGM includes many subtle effects, with degen- 
eracies arising from fluctuations in density, temperature, 
and ionization fraction, as well as spectral hardness of 
the ionizing sources. Previous estimates of z T « 3 were 
based on indirect measurements of temperature, ioniza- 
tion, and average optical depth of Lya lines of H I, and on 
apparent shifts in Si IV/C IV absorption-line ratios aris- 
ing from changes in the 4-ryd continuum radiation field 
as He II was ionized to He III. Early estimates of the 
IGM temperature evolution and equation of state from 
H I line widths suggested a change at z « 3 (Ricotti, 
Gnedin, & Shull 2000; Schaye et al. 2000). Theuns et al. 
(2002) analyzed the average optical depth of (H I) Lya 
lines in 1061 quasar spectra, using low-resolution data 
from the Sloan Digital Sky Survey (Bernardi et al. 2003). 
From a weak dip in t c s(z) at z « 3.1, these authors sug- 
gested that He II reionization started at redshift z ~ 3.4 
and lasted until z « 3.0. Several recent studies with 
higher-resolution spectroscopic data have not confirmed 
the r eff feature (McDonald et al. 2005; Kim et al. 2007). 
Faucher-Giguere et al. (2008) found a feature in r c ff (H I) 
near z = 3.2, in a sample of 86 high-resolution, high- 
S/N quasar spectra taken with Keck (ESI, HIRES) and 
Magellan (MIKE). However, they question the reioniza- 
tion explanation of the depth and narrow redshift extent. 
They also note that at least three physical effects could 
produce a feature in in t c s (z): changes in IGM tempera- 
ture, electron number density, and ionizing background. 
Recent numerical simulations (Bolton et al. 2009; Mc- 
Quinn ct al. 2009; Lidz et al. 2010) also question whether 
He II reionization is responsible for the proposed feature 
in t c s(z) at z « 3.2 ± 0.2. 

The most direct estimates of He II reionization come 
from far-ultraviolet spectra of He II Lyman-a absorption 
toward high-redshift AGN. Because of strong Galactic 
interstellar absorption in the H I Lyman continuum, the 
He II A304 line is only observable in the far ultravio- 
let, at redshifts z > 2 when the line shifts longward of 
912 A. A small number of bright AGN at z = 2.9-3.3 have 
been observed with the Hubble Space Telescope (HST) 
and the Far Ultraviolet Spectroscopic Explorer (FUSE). 
These data show moderate optical depths in the He II 
line, ranging up to th c ii = 3 — 4, with a gradual re- 
covery of transmission (rneii < 1) at z < 2.7. The 
best He II absorption data have been acquired toward 
three high-redshift AGN: Q0302-003 (Jakobsen et al. 
1994; Hogan, Anderson, & Rugers 1997; Heap et al. 
2000); and HE 2347-4342 (Reimers et al. 1997; Kriss 
et al. 2001; Smctte ct al. 2002; Shull et al. 2004; Zheng 
et al. 2004); and HS 1700+6416 (Davidsen, Kriss, & 
Zheng 1996; Fechner et al. 2006). From an analysis of 



the He II optical-depth evolution in moderate-resolution 
FUSE data toward HE 2347-4342, Shull et al. (2004) 
suggested that the reionization epoch of He + occurred 
at z r = 2.8 ± 0.2. 

Observations of He II with HST and FUSE also con- 
strain the thermodynamic state of the IGM and the spec- 
tra of ionizing sources. For example, the high abun- 
dance ratio, N(He II)/N(H I) w 50-100, observed toward 
HE 2347-4342 and HS 1700+6416 (Kriss et al. 2001; 
Smette et al. 2002; Shull et al. 2004; Zheng et al. 2004; 
Fechner et al. 2006) requires a soft ionizing radiation 
field, F v oc ^~ Qs , with mean spectral index (a s ) ~ 1.8 
(Fardal et al. 1998). As observed in the QSO rest-frame, 
the EUV (1-2 ryd) spectral indices of quasars exhibit a 
wide range (Telfer et al. 2002; Scott et al. 2004). More 
intriguing are the spectroscopic observations of He II and 
H I taken by FUSE, the Keck Telescope, and the Very 
Large Telescope (VLT). Kriss et al. (2001) and Shull et al. 
(2004) report variations in r] = 20 - 200 over 2-10 Mpc 
scales (Az « 0.002 — 0.01) which imply spatial fluctu- 
ations in the radiation field at the ionization edges of 
H I (1 ryd) and He II (4 ryd). Most of these variations 
probably arise in the 4 ryd continuum, whose photons 
have greater penetrating power. These flux differences 
at 1 ryd and 4 ryd are then amplified by absorption and 
reprocessing by the IGM (Haardt & Madau 1996; Fardal 
et al. 1998). 

In this paper, we present new, high-quality, far-UV 
spectroscopic observations of intergalactic He II absorp- 
tion toward HE 2347—4342, using the Cosmic Origins 
Spectrograph (COS) on HST (Green et al. 2010; Oster- 
man et al. 2010). This target is of considerable inter- 
est as a bright (V = 16.1, z ~ 2.9) continuum source, 
previously observed by the FUSE, Keck, and VLT spec- 
trographs to study intergalactic Lya absorption in both 
H I and He II. HE 2347-4342 is the first of three AGN 
targets scheduled for COS guarantecd-time observations 
of the He II reionization epoch. The data were ac- 
quired with both the moderate-resolution G130M grating 
(R w 18, 000, A = 1135 - 1440 A) and the low-resolution 
G140L grating (R « 1500, A = 1030 - 2000 A). In the 
303.78 A line of He II, these wavelength bands allow us 
to probe redshifts down to ZHcii = 2.735 (G130M) and 
ZHeii ~ 2.39 (G140L). The COS data enable many im- 
provements compared to previous studies. First, the high 
far-UV throughput of the COS/G130M grating provides 
higher signal-to-noise (S/N) and better photometric ac- 
curacy. Second, the low background of the COS detectors 
allows us to detect very weak flux-transmission through 
the IGM, characterize low flux levels in He II absorp- 
tion troughs, and probe regions of high optical depth, 
THcii 5. Finally, low-resolution G140L spectra (1000- 
2000 A) allow us to study the recovery of the He II op- 
tical depth at z < 2.7 and characterize the continuum 
longward of the He II edge. 

In §2 we discuss the observations and data reduction 
techniques for both G130M and G140L gratings. In §3 
we display the fitted quasar continuum, Fx oc A _3 0±01 , 
which we use to extrapolate below the He II edge at 
1186.2 A. From this continuum and the transmitted 
fluxes, we derive He II optical depths in the Lya forest 
and analyze the fluctuating absorption. We observe min- 
imal effects of local photoionization expected from prox- 
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imity to this very luminous QSO. In the He II Lya forest 
at 2.4 < z < 2.9, we see narrow (2-5 A, 12-30 Mpc of 
comoving radial distance) windows of partial flux trans- 
mission, which we interpret as "patchy He II reioniza- 
tion", interspersed with broad He II absorption troughs 
(4-10 A intervals, 25-60 Mpc) with high optical depths, 
THcii > 5. At A < 1130 A (zhoII < 2.72) the mean opac- 
ity drops to (thcIi) ^ 2, consistent with the overlap of 
the He II I- fronts. Because He II is much more abundant 
than H I, its opacity through the cosmic web remains 
high. 

In §4 we summarize our observations and their impli- 
cations for the He II reionization epoch. With reason- 
able assumptions about the physical state of the IGM at 
z = 2 — 3, borne out by studies of H I and He II Lya- 
forest absorption and cosmological simulations, we can 
probe the He II reionization epoch down to fractional 
ionization ZHcii > (0.01)(0.1/5jj)(rH e ii/5) in underdcnse 
regions (Sh = 0.1 — 1) relative to the mean baryon density 
of the cosmic web. This ionization level provides much 
greater sensitivity to conditions at z ps 2.8 than possible 
with H I at z ps 6 (Fan et al. 2006). Coupled with galaxy 
surveys, cosmological simulations, and photoionization 
modeling, these UV spectra will help us understand the 
propagation and overlap of the helium ionization fronts. 

2. HST-COS OBSERVATIONS AND DATA 
REDUCTION 

2.1. Wavelength Scale and Fluxes 

HE 2347-4342 was observed with the short- 
wavelength, medium-resolution far-UV mode of 
HST-COS (G130M) on 2009 November 05 for a 
total of 28,459 s. Descriptions of the COS instrument 
and on-orbit performance characteristics are in prepa- 
ration (Green et al. 2010; Osterman et al. 2010). In 
order to achieve continuous spectral coverage across 
the G130M bandpass (1135 A < A < 1440 A) and to 
minimize fixed-pattern noise, we made observations at 
two central wavelength settings (1291 A and 1300 A) 
with four focal-plane offset locations in each grating 
setting (FP-POS = 1, 2, 3, 4). This combination of 
grating settings ensures the highest signal-to-noise 
observations at the shortest wavelengths available to the 
G130M mode (zueii ~ 2.73 — 3.00) at a resolving power 
of R = (A/AA) w 18, 000 (Ad FW hm « 17 km s" 1 ). 

In addition, COS observed HE 2347-4342 in the low- 
resolution G140L mode for 11,558 s on 2009 November 
06. The G140L Segment-A data provide some over- 
lap with the G130M observations between 1240-1440 A, 
and they extend our low-resolution spectral coverage 
(R w 1500) out to A 2000 A. From G140L segment-B 
observations and custom data-processing steps described 
in Section 2.2, we extend our reliable spectral coverage 
over part of the G140L bandpass, between 1030-1175 A. 
Table 1 lists the relevant COS observational parameters 
for the data sets used in this work. 

All observations were centered on the QSO tar- 
get, HE 2347-4342 at J2000 coordinates (R.A. = 
23 h 50 m 34.23 s , Deck = -43°25'59.80"). COS per- 
formed NUV imaging target acquisitions with the MIR- 
RORA/PSA mode, yielding good centering in the COS 
aperture to maximize throughput and resolving power. 



TABLE 1 
HE 2347-4342 Observing Log 



Date 


COS Mode 


Ac (A) 


FP-POS 


Exp 




2009 Nov 05 


G130M 


1291 A 


1,2,3,4 


9 


16535 


2009 Nov 05 


G130M 


1300 A 


1,2,3,4 


9 


11924 


2009 Nov 06 


G140L 


1230 A 


1,2,3,4 


i 


11558 



The G130M data were processed with the COS calibra- 
tion pipeline, CALCOSQ v2.11b, and combined with a 
custom IDL coaddition procedure. The analog nature of 
the COS micro-channel plate makes it susceptible to tem- 
perature changes. Electronically injected pulses (stims) 
at opposite corners of the detector allow for the track- 
ing and correction of any drift (wavelength zero-point) 
and/or stretch (dispersion solution) of the recorded data 
location as a function of temperature. At the begin- 
ning of the observations of HE 2347—4342, the detector 
proved to be substantially colder than the nominal op- 
erating temperature of the micro-channel plate. This 
allowed one of the electronic stim-pulses to accumulate 
at the bottom of the detector, compromising the wave- 
length solution for the first four G130M exposures. The 
detector warmed to the normal operating range for the 
remaining G130M observations. We reprocessed the ear- 
lier exposures using extrapolated stim positions to ensure 
a consistent wavelength solution over the course of the 
observations. The G140L data were unaffected. 

Flat-fielding, alignment, and coaddition of the pro- 
cessed COS/G130M exposures were carried out using 
IDL routines developed by the COS GTO team specif- 
ically for COS FUV dat£0- First, the data were cor- 
rected for fixed-pattern instrumental features. Although 
attempts at a true "flat-fielding" of COS data show 
promise, the technique is not yet sufficiently robust to 
improve data of moderate S/N. The quantum efficiency 
is improved by the presence of a series of wires, called 
the quantum-efficiency grid, placed above the detector. 
These wires create shadows in the spectrum that are 
removed during data reduction. An ion-repellor grid 
reduces the background rate by preventing low-energy 
thermal ions from entering the open-faced detector. We 
are able to correct the narrow shadows of ^15% opac- 
ity arising from the grid wires. A one-dimensional map 
of grid-wire opacity for each detector was shifted from 
detector coordinates into wavelength space and divided 
from the flux and error vectors. Exposure time in the lo- 
cation of grid wires was decreased by ~ 0.7, giving these 
pixels less weight in the final coaddition. We also modi- 
fied the error vector and local exposure time at the edges 
of the detector segments to underweight the flux contri- 
butions from these regions. With a total of eight different 
wavelength settings in the final co-addition, any residual 
instrumental artifacts from grid-wire shadows and detec- 
tor segment boundaries should have little effect on the 
final spectrum. 

See the HST Cycle 18 COS Instrument Handbook for more 

details: 

http : //www. stsci . edu/hst/cos/documents/handbooks/current/cos 
4 See http://casa.colorado.edu/~danforth/costools.html 

for our coaddition and ilat-rieiding aigorithm and additional dis- 
cussion. 
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We performed two sets of wavelength-scale shifts. 
First, the G130M exposures were aligned with each other 
and interpolated onto a common wavelength scale. Then, 
we transferred the G130M wavelengths to a ground-based 
reference frame of the H I data, taken by D. Reimers 
with the VLT/UVES (Ultraviolet Visual Echelle Spec- 
trograph) in ESO program 071.A-0066. These data were 
used in the He II analysis paper of FUSE data (Shull 
et al. (2004). One exposure in each COS detector was 
picked as a wavelength reference, and the remaining ex- 
posures were cross-correlated with it. The wavelength 
region of cross-correlation for each case was picked to in- 
clude a strong interstellar absorption feature, although 
shifts were typically a resolution element or less. Resid- 
ual errors on the target centering within the COS aper- 
ture caused a systematic shift between the strong inter- 
stellar lines and their expected LSR velocities. The COS 
wavelength solution is relatively accurate, but it may suf- 
fer from linear offsets from an absolute frame. Since our 
later analysis involves comparing the rapidly- varying H I 
(optical) and He II (UV) profiles, wavelength alignment 
is critical. The only absorption features in the normal- 
ized VLT data clearly not associated with the Lya forest 
are a pair of narrow absorbers at 4787.12 A and 4795.08 
A, consistent in relative strength and separation with the 
C IV doublet at 1548.195 and 1550.770 A, redshifted to 
z = 2.0921. At this redshift, COS/G130M covers 384- 
466 A. The C III A386.20 line at 1194.1 A in the COS data 
corresponds to the C IV doublet in the optical data, with 
similar line profiles. We therefore shift the COS data by 
+0.07 A (approximately one G130M resolution element) 
to align with the optical spectra. 

The flux at each position was taken to be the exposure- 
time weighted mean of flux. Since exposure time was re- 
duced in certain wavelength locations, as noted above, 
pixels near detector edges and where grid-wire shad- 
ows were removed received less weight than those in 
less suspect locations. The combined G130M data show 
S/N ~ 25 per 7-pixel resolution element, sufficient to 
detect narrow absorption features near 1200 A down to 
equivalent widths W\ > 10 mA at 4er significance. The 
G140L observations were originally processed with CAL- 
COS 2.11b. However, the A < 1200 A data from this 
processing is not appropriate for scientific analysis. We 
describe our custom processing in Section 2.2. 

Figure 1 shows the G130M/G140L continuum from 
1200 A to 2000 A, which we have fitted and extrap- 
olated below the He II edge at 1186.2 A. The G140L 
wavelengths were shifted to match G130M, with no cor- 
rections needed to match G130M and G140L fluxes. 
These curves were de-reddened using color excess, E(B — 
V) = 0.022, corresponding to column density log Am = 
20.06 derived by fitting interstellar Lya absorption. We 
adopted a standard interstellar ratio, Nhi = (5.3 x 
10 21 cm" 2 mag- 1 ) J B( J B - V), from Shull & Van Steen- 
berg (1985) and Diplas & Savage (1994) and used the UV 
selective extinction curve of Fitzpatrick & Massa (1999). 
Our derived extinction, E(B — V) = 0.022, lies inter- 
mediate between previous values of E(B — V) = 0.0387 
(Reimers et al. 1997) and E(B - V) = 0.014 (Smette 
et al. 2002). Figure 2 shows the G130M spectrum 
(A < 1208 A) with prominent He II absorption troughs, 



TABLE 2 

HE 2347-4342 Sight Line He II Absorption 



A.z 



A c (A) 



r HoII 



COS Mode 



Notcs b 



2.90151 
2.89765 
2.88721 
2.87445 
2.84700 
2.77471 
2.71324 
2.69082 
2.66414 
2.65061 
2.63785 
2.61814 
2.59804 
2.57059 
2.55000 
2.52999 
2.49327 
2.47001 
2.43683 



0.00400 
0.00800 
0.00823 
0.01317 
0.03456 
0.05267 
0.00823 
0.01646 
0.00494 
0.00494 
0.01317 
0.00988 
0.00988 
0.01317 
0.02008 
0.02008 
0.01975 
0.02008 
0.01811 



1185.20 
1184.03 
1180.86 
1176.98 
1168.64 
1146.68 
1128.01 
1121.20 
1113.09 
1108.98 
1105.11 
1099.12 
1093.01 
1084.67 
1078.42 
1072.34 
1061.19 
1054.12 
1044.04 



3.3 ±0.25 
3.0 ± 0.50 
1.8 ±0.25 
1.5 ±0.25 
2.3 ±0.25 
1.3 ±0.50 
2.0 ±1.0 
3.5 ± 1.0 
1.0 ±0.25 
2.5 ±0.50 
2.8 ± 1.0 
1.5 ±0.50 
1.7 ± 0.50 



G130M 
G130M 
G130M 
G130M 
G130M 
G130M 
G140L 
G140L 
G140L 
G140L 
G140L 
G140L 
G140L 
G140L 
G140L 
G140L 
G140L 
G140L 
G140L 



Z04-3A 
Z04-3CDE 
Z04-3HIJ 



a Average Hell optical depth THell was computed from flux trans- 
mission averaged over redshift bins z c ± Az at centroid wavelengths 
A c . Errors in THell are conservatively rounded to nearest 0.25. 
b These three He II complexes noted Z04-3n were identified in Table 
3 of Zheng et al. (2004). 



numerous flux-transmission windows shortward of the 
He II edge, and the two proposed QSO systemic red- 
shifts, z Q so = 2.885 ± 0.005 (Reimers et al. 1997) and 
zqso — 2.904±0.002 (this paper). This issue is discussed 
further in Section 3.1. 

The strong He II features are listed in Table 2, with 
the He II optical depths calculated from the depth of flux 
transmission, th c ii = h\(Fo/F\), relative to the extrap- 
olated continuum F . Uncertainties in rneii were com- 
puted as root-mcan-square averages over redshift bins, 
Az, ranging from 0.004 to 0.02 at wavelength centroids 
A c . We conservatively rounded errors to the nearest 0.25 
in Tfjeii, using uncertainties arising from converting er- 
rors in flux transmission to errors in optical depth. In 
regions of lower S/N in the G140L grating (z < 2.6), we 
smoothed over two resolution elements (~ 2 A) to re- 
duce the noise before the optical depth and error were 
evaluated. The G140L-B spectra shown in Figure 1 were 
smoothed over three resolution elements. 

2.2. Data Processing: COS G140L Segment B 

During the Servicing Mission-4 Observatory Verifi- 
cation (SMOV-4) period, it was discovered that the 
MgF2/Al mirrors of H ST have retained approximately 
80% of their prc-flight reflectivity at A < 1175 A (Mc- 
Candliss et al. 2010). This result has opened the door 
for use of the short-wavelength response of COS, in par- 
ticular the G140L mode, to perform spectroscopic obser- 
vations at wavelengths inaccessible to previous HST in- 
struments (A < 1130 A). Sec Appendix A of France et al. 
(2010) for preliminary work with this mode. The CAL- 
COS pipeline processing of G140L segment B (400 A 
< A < 1175 A) is not sufficiently mature to produce one- 
dimensional spectra appropriate for scientific analysis. 

We therefore processed the data with a newer version 
of the CALCOS pipeline (v2.12), employing specially 
created reference files for flux and wavelength calibra- 
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1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 

Wavelength (A) 

Fig. 1. — Overview of the COS/G130M and G140L data. Data are smoothed by 3 resolution elements and shown as observed (uncorrected) 
flux. The de-reddened data were fitted to a power law, Fx = F (A/A )- as , with a a = 3.0±0.1. Error bars below A = 1215 A reflect ±2% 
uncertainty in continuum level and ±0.1 uncertainty in ct s . This continuum is extrapolated below the Hell edge (1186 A) to derive Hell 
optical depths. Red curve shows continuum fit of Kriss et al. (2001) with ct s = 2.4. Our steeper fit arises primarily from the higher quality 
COS data (1200-2000 A) and, to lesser extent, from our larger adopted extinction, E(B — V) = 0.022 rather than 0.014, extrapolated to 
the far-UV. 



tions in the short-wavelength segment-B setting. The 
custom reference files that were created to supplement 
the CALCOS 2.12 release were a first-order dispersion 
solution for A < 1150 A and a flux calibration curve. 
The wavelength calibration was created by fitting line 
centers to the pressure-broadened Lyman series lines ob- 
served in COS calibration spectra of white dwarfs (Mc- 
Candliss et al. 2010). While G140L segment A uses a 
second-order polynomial wavelength solution, the rela- 
tively small amount of physical detector area used by 
G140L segment B allows a reasonable fit with a linear so- 
lution. We find a dispersion of AA = —79.8 mA pixel -1 , 
and a long- wavelength zero-point of 1245.26 A. To align 
with G130M data, we shifted G140L-B wavelengths by 
+3.9 A. For example, 1245.26 A (old G140L) becomes 
1249.16 A (new G140L). The G140L-A wavelengths were 
shifted by only —0.4 A (~ 0.5 resolution element). 

The calibration curve, translating count rate 
(counts s _1 ) into flux (erg cm -2 s _1 A -1 ) was 
created using the G140L segment-B effective area curve 
derived by McCandliss et al. (2010). We caution that 
there may still be significant uncertainty in this flux 
calibration. At A < 950 A, we should be within ~ 50% 



of the correct value. The accuracy is likely —20-25% 
from 1025-1110 A and -10-15% from 1120 A longward. 
Continued testing of G140L-B observations will refine 
this calibration. This custom reduction produces repro- 
cessed one- and two-dimensional spectral data products, 
which can be coadded following the procedure described 
above and in Danforth et al. (2010). The coadded 
data were used to create a low-resolution (R » 10 3 , 
AA - 1 A) spectrum of HE 2347-4342 at A < 1174 A. 
In principle, the COS/G140L wavelength coverage 
extends to very low wavelengths (~ 400 A). However, 
the combination of low instrumental effective area, 
A eff ~ 10 cm 2 at A < 1030 A (McCandliss et al. 2010) 
and interstellar absorption by atomic and molecular 
hydrogen limit the working wavelength coverage. For 
the He II absorption toward HE 2347—4342, we present 
G140L data from 1030 A < A < 1174 A, corresponding 
to zhoII ~ 2.39 - 2.87. We cross-calibrated the G140L 
segment-B spectrum with the short-wavelength data 
from G130M, finding that the flux calibration is better 
than 10% in the overlap region. The data were then 
used to extend our investigation of the He II ionization 
state in the IGM down to z Hc ii = 2.39 (1030 A). 
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Fig. 2. — Detail of the G130M data showing Hell absorption and transmission windows. The two proposed QSO systemic redshifts 
(zqso = 2.885 and 2.904) and our extrapolated continuum are marked. Milky Way interstellar lines of NI, Sill, Si III appear longward 
of 1190 A, and strong Hell absorption is seen shortward of 1186.26 A. We use the CIII (z = 2.0921) IGM absorption line to align the 
COS and VLT spectra. Windows of flux transmission appear at 1138, 1144, 1154, 1160, 1174, 1179, 1183, and 1184 A. At 1174 A we see 
essentially 100% transmission, and yet we also observe three troughs of strong Hell absorption (Tu e ii > 5) over large redshift intervals 
(and comoving radial distances) between * = 2.751 - 2.807 (61 Mpc), z = 2.823 - 2.860 (39 Mpc), and z = 2.868 - 2.892 (25 Mpc). 



Figure 3 provides a comparison among three spectro- 
graphs: COS/G130M, COS/G140L, and FUSE (LiFl 
and LiF2 channels). There is a fair amount of agree- 
ment in fluxes among these instruments, although the 
FUSE backgrounds at 1160-1187 A are obviously incor- 
rect. This comparison gives us confidence in using G140L 
data down to 1030 A. 

3. RESULTS AND ANALYSIS 

3.1. Systemic Redshift of QSO 

The systemic redshift of HE 2347-4342 is of inter- 
est, both to fix the expected location of the He II edge 
and to infer the metagalactic radiation field from the 
QSO proximity effect (Dall'Aglio, Wisotzki, & Worseck 
2008a, b). A significant puzzle for this sight line is the 
6 A offset between the extent of strong He II absorp- 
tion (out to 1186 A) and the expected He II edge (1180 
A) based on an uncertain QSO emission-line redshift. 
Reimers et al. (1997) found z cm = 2.885 ± 0.005, based 
on the weak O I A1302 emission line, while Scannapieco 
et al. (2006) quoted z cm = 2.871 with no error bar. 
Emission-line redshifts were also measured with UVES at 
the VLT/UT2 by Dall'Aglio et al. (2008b), who quoted 
z ero = 2.886 ± 0.003 based on the O I + Si II blend and 
assuming a rest wavelength A rost = 1305.77 for the blend. 
Given the solar abundance ratio (O/Si)© » 15, this blend 
is probably dominated by three lines of O I, with a statis- 



tically weighted centroid at 1303.49 A (Morton 2003). In 
total, the blend consists of three O I lines decaying to the 
ground-state fine-structure levels of O I [(2p 3 3s) 3 Si — > 
(2p 4 ) 3 P 2 ,i, ] at 1302.17, 1304.86, 1306.03 A and two 
lines of Si II [(3s 3p 2 ) 2 S 1/2 ->• (3s 2 3p) 2 P lA3/2 ] at 

1304.37 and 1309.28 A. Dall'Aglio et al. (2008a) mea- 
sured the redshifts of the strongest emission lines to be 
z ein = 2.877±0.003 (Lya) and 2.861±0.003 (C IV). How- 
ever, there are well-known shifts in emission-line redshifts 
between the QSO systemic redshift and the "higher ion- 
ization lines" (e.g., Lya, C III A977, C IV A1549). For 
instance, 1000-3000 km s _1 offsets are commonly seen 
in the emission-line spectra of AGN (Espey et al. 1989; 
Corbin 1990). These issues provide systematic uncertain- 
ties in determining z cm , either from the weak O I/Si II 
blend or from higher ionization lines. Thus, the precise 
QSO systemic redshift remains uncertain. 

Because of these problems with emission-line redshifts, 
we propose to use IGM absorption to determine the QSO 
redshift. We suggest that the QSO systemic redshift is 
zqso = 2.904 ± 0.002 (Figure 4) based on the extent of 
the H I and He II absorption. This redshift is chosen as 
the midpoint between the centroid of the strong absorp- 
tion feature seen in both H I and He II and the red edge 
of He II absorption. Figure 4 provides an overview of 
the COS/G130M spectrum from 1100-1190 A, together 
with a blowup of the region from 1173 A to the He II 
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Fig. 3.— Comparison of the COS/G130M (black), COS/G140L (blue), and coadded FUSE/LiFl+LiF2 data (green). The COS/G130M 
data are binned by 3 pixels (half a resolution element). The FUSE data have been binned to one resolution element, and COS/G140L 
data are presented as observed (unbinned). The FUSE background solutions appear to be adequate at 1110-1150 A, but poor between 
1160-1187 A. The flux and wavelength solutions for the two instruments are reasonably consistent, providing confidence in using G140L 
data at A > 1030 A. 



edge near 1186 A. We see absorption extending well be- 
yond the wavelengths (1180.2-1180.5 A) corresponding 
to the QSO emission-line redshifts, z — 2.885 ± 0.005 
and 2.886±0.003, suggested by Reimers et al. (1997) and 
Dall'Aglio et al. (2008b) respectively. The absorption be- 
tween 1181-1186 A suggests that the QSO actually lies 
at z w 2.90. Interestingly, we see no ionization effects 
from proximity to the QSO, a topic discussed further in 
Section 3.6. 

The He II absorption profile shows an absorption fea- 
ture at 1186.5 A, with optical depth r w 0.7, equivalent 
width 95 mA, and FWHM = 0.14 A. This may be proxi- 
mate He II absorption or absorption from an intervening 
system. The spectrum near 1186.5 A shows a clear flux 
minimum and recovery toward shorter wavelengths typ- 
ical of weak absorption lines. A distinct He II absorber 



would appear as a broader "shelf on the side of the 
THcii > 1 trough. We do not make a specific line identifi- 
cation for the feature. It does not seem to be low-redshift 
Ly/3 nor any of the usual IGM metal-line systems (Dan- 
forth & Shull 2008). 

Reimers et al. (1997) attribute the H I and He II 
absorption at z > 2.885 to a multicomponent system 
of associated absorbers between 1181-1186 A (redshifts 
z = 2.891 to z = 2.904). Many of these systems have 
strong metal lines of C IV, N V, and O VI, and a few 
have anomalous ratios of He II and H I Lya absorption. 
These observations suggest that the gas might be exposed 
to high fluxes of ionizing radiation and affected by nucle- 
osynthetic anomalies in nuclear outflows. Smette et al. 
(2002) and Fechner et al. (2004) both noted evidence for 
"line-locking" in C IV and N V (although the latter group 
found no statistical evidence compared to simulated line 
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Fig. 4.— (Top Panel) Normalized COS/G130M and G140L flux (black) and VLT flux (red) vs. redshift (z = 2.4 - 2.9, top legend) and 
wavelength (bottom legend). HI wavelengths are divided by 4 for comparison with Hell, and two proposed QSO rcdshifts are noted by 
filled black circles: z cm = 2.885 ± 0.005 (Reimers et al. 1997) and z = 2.904 ± 0.002 (this paper). Lower two panels show VLT/UVES 
data (R = 44000) from D. Reimers, convolved with COS line spread function to match COS resolution: R = 18,000 (A > 1134 A) and 
R = 1800 (A < 1134 A). Unconvolved HI data are shown by dotted red curves. Troughs of strong Hell absorption (1175.2-1179.1 A, 
1179.5-1182.5 A) are separated by small windows of flux transmission. Features at 1181-1186 A are associated absorbers (Reimers et al. 
1997), but with insufficient column density to block AGN ionizing flux. 
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lists). Line- locking occurs when two absorbers are sep- 
arated by a velocity equal to the separation of the two 
doublet lines, resulting from radiation pressure on ma- 
terial physically near the QSO (Weymann, Carswell, & 
Smith 1981). These radiative forces will drive gas out- 
ward from the QSO, whereas the redshifted associated 
absorbers would have infall velocity ~ 1500 km s -1 rel- 
ative to the QSO if zqso = 2.885. Thus, our proposed 
rcdshift, zqso = 2.904, provides a more plausible expla- 
nation. 

The uncertainty in QSO systemic redshift illustrates 
the need for infrared observations. Unfortunately, the 
redshifted hydrogen Balmer lines fall in difficult spectral 
bands (Ha at 2.563 /im, H/3 at 1.898 /um), as does the 
[O II] A3727 doublet at 1.455 and 1.456 /zm. The forbid- 
den lines of [O III] A5007, 4959 occur at more promising 
wavelengths, 1.955 /xm and 1.936 /an. Thus, infrared 
spectroscopic searches for the [O III] forbidden lines may 
offer the best opportunity to resolve the issue of zqso 

3.2. Expected He II Gunn-Peterson Absorption 

The standard GP optical-depth formula (Fan et al. 2006) 
for species (i = H I or He II) is 



are: 



n(z) = ( ^ 

TO„C 



K fi n i 



(1) 



where (ire 2 /m e c) = 2.654 x 10~ 2 cm 2 Hz is the integrated 
classical oscillator cross section, /, = 0.4162 is the ab- 
sorption oscillator strength for Lya transitions of H I and 
He II at rest wavelengths A» = 1215.67 A and 303.78 A, 
and rii is their number density (cm -3 ). The Hubble ex- 
pansion parameter is H(z) = H [fl m (l + z) 3 + ^a] 1//2 in 
a flat Friedman cosmology (fi m + fi"A = 1)- All cosmolog- 
ical parameters and the resulting coefficients are based 
on WMAP-7 concordance values of CMB modeling from 
Komatsu et al. (2010), n m h 2 = 0.1349 ± 0.0036, ft A = 
0.728±°;°^, and Vt b h 2 = 0.02260 ± 0.00053, and Hubble 

constant H = (100 km s _1 Mpc^ 1 )h with h = 0.70. To 
set the ion densities rii, we define oihi and XHeii as the 
abundance fractions of H I and He II, respectively. (The 
quantity xhi is often called the neutral fraction.) The 
factors Si are the overdensities of H (at z « 6) and He 
(at z « 3) relative to the mean co-moving baryon den- 
sity, (p b ) = (fifePcr) = 4.25 x 10~ 31 g cm" 3 . With a pri- 
mordial helium abundance Y = 0.2477 ± 0.0029 by mass 
(Peimbert et al. 2007) and y = 074) /(l - Y) = 0.0823 
by number, this baryon density corresponds to H I and 
He II number densities, 



n H i = (6.55 x 10 5 cm 3 ) x m Sh 



(1 + *) 



n H eii = (l-01 x 10 6 cm 3 ) a; He ii fee 



(1 + *) 



(2) 



.(3) 



In the high-z limit, O m (l + z) 3 3> ^a, we can express 
the Hubble parameter as H(z) « H flm 2 (l + z) 3 / 2 « 
(294 km s" 1 Mpc _1 )[(l + z)/4)] 3 / 2 scaled to the epoch 
z w 3 relevant for the beginning of He II reionization. 

Thus, the Gunn-Peterson optical depths for H I (scaled 
to the z = 6 epoch) and He II (scaled to the z = 3 epoch) 



Thi(2) = (3.99 x 10 5 ) 
t Ho ii(z) = (3.55x 10 3 ) 




xm fe , (4) 



XHell feell ■ (5) 



Very small abundance fractions, xhi ~ 10 -4 and Xfjoii ~ 
10~ 2 , will saturate the Gunn-Peterson Lya absorption 
in H I and He II, with optical depths Tj 3> 10, even in 
under-dense regions with Si = 0.1. It is primarily for this 
reason that the hydrogen reionization epoch is so diffi- 
cult to measure in Lya absorption at z > 6 (Fan et al. 
2006). However, He II absorption at z < 3 can probe 
50-100 times higher abundance fractions in XHcii- For 
the maximum observable values of He II optical depth 
(thcII = 5) and with typical densities in low-density fil- 
aments and voids, the intervals of high He II absorption 
at z = 2.8 ±0.1 correspond to He II fractions, 



ZHell = (0-015) 



THoII 



(1+f) 

3.8 



-3/2 



Skc 

0.1 



(6) 



These overdensities, <5h — 0.1, arc appropriate for the 
low-column-density (log A^hi _ 13.5) regions of the Lya 
forest. 

Figure 5 shows the He II optical depths for 2.4 < z < 
2.9, using COS data from both G140L and G130M and 
FUSE. To construct this figure, we reconstructed "hor- 
izontal error bars" representing the redshift bins (Az) 
used in these papers. We used unpublished Az and un- 
certainties in optical depths THeii from Shull et al. (2004) . 
From Zheng et al. (2004), we used values of Az and THeii 
from their Figure 3, with 25% uncertainties in t Hell* We 
find reasonable agreement with the FUSE optical depths 
(Shull et al. 2004) down to z = 2.4, although some of the 
G140L fluxes are higher than those in the FUSE exten- 
sion to z < 2.4 (Zheng et al. 2004). The smaller FUSE 
optical depths may arise from uncertain background sub- 
tractions. Defining the redshift ranges for these data 
points is difficult for optical depths that vary rapidly over 
AzwO.01. 

Future COS spectra of fainter He II-GP targets (of- 
ten 19 th magnitude AGN) will have low S/N ratios and 
will probably not show the fluctuations in optical depths 
found here. At best, they will measure broad-band op- 
tical depths, (th c ii). For HE 2347-4342, we list below 
the He II and H I optical depths averaged over Az = 0.1 
windows, together with the standard deviations of their 
distributions: 

. (z = 2.4-2.5) (thi) = 0.12±g;g| , (r He ii) = l.llto.49 
. (z = 2.5-2.6) (thi) = 0.12™ , (t Ho ii) = 1.18tg;& 

. (z = 2.6-2.7) (thi) - 0.26™ , (t Hc ii) - l-37t^ 3 

There may be a slight trend to higher r\ at lower redshifts, 
as we discuss in Section 3.4. Considering the ratios of the 
broad-band averages, we find that (77) rises from 21± 33 , 
at z = 2.6 - 2.7 to values 39±|t at z = 2.5 - 2.6 and 



37; 



-72 
-34 



at z = 2.4 — 2.5. Despite the large spreads of 
the optical-depth distributions, the means are fairly well 
determined. 
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Fig. 5. — Absorption optical depths, th ii, vs. redshift (bins listed in Table 2) for this paper (black points: COS G130M and G140L) and 
from previous FUSE measurements (Zheng et al. 2004; Shull et al. 2004). COS optical depths are generally higher than those found by 
FUSE. Theoretically predicted (1 + z) 3,5 redshift scaling (Fardal et al. 1998) is shown as a dotted line. High COS optical depths, thcII > 5, 
are plotted as lower limits, based on zero-level flux accuracy implied by RMS variations, ap = 2.7 X 10~ 17 ergs cm" 2 s _1 A -1 , at the 
bottom of troughs at 1145-1155 A and 1162-1173 A. 



Figure 6 shows examples of the He II and H I ionization 
structures found in our cosmological simulations (Smith 
et al. 2010), computed with Enzo, an N-body, hydro- 
dynamical code modified to compute time-dependent H 
and He ionization. One can see the projected topology 
of the cosmic web, the baryon overdensity, the ioniza- 
tion fractions of H I and He II at z = 2 — 3, and their 
ratio, Xneii/xm- Observationally, we measure the ra- 
tio of column densities, r\ = N(He II)/N(H I), which we 
approximate as rj s=s 4th c ii/'?hi- The factor of 4 arises 
from the larger He II frequency bandwidth, with the 
absorption-line cross section, o\ cx /A, scaling as oscil- 
lator strength times wavelength. Our IGM simulations 
at z « 2.5 show a mean value of the ionization ratio, 
^Heii/^Hi ~ 500, corresponding to (rj) « 40 for primor- 
dial He/H = 0.0823 by number. We have confirmed that 
a substantial amount of observable He II, with XHcii = 
0.001-0.01, comes from ionized gas with logT = 4.0 — 5.3 
(much of it at logT = 4.3) and logfe = -1 to +2. This 
He II gas is primarily confined to photoionized filaments 
in the cosmic web, with some contribution from colli- 
sional ionization (at T > 10 5 K). Further analysis of 



these simulations will appear in a later paper. 

3.3. He II and H I Optical-Depth Ratios 

The H I and He II Lya-forest absorbers have tradi- 
tionally been assumed to arise in photoionized gas, irra- 
diated by metagalactic ionizing (EUV) radiation of spe- 
cific intensity J„ = Ji(v / i^)~ Qs , where i refers to H I 
(hv l = 13.6 eV) or He II {hv l = 54.4 eV). For these 
hydrogenic species, whose photoionization cross sections 
above threshold hvi scale as er„ = <7i(i//i/i) , the pho- 
toionization rates = (incri/h) Ji/[oti +3). The param- 
eters Oii are local spectral indices of the ionizing back- 
ground at 1 ryd and 4 ryd, respectively, and provide mi- 
nor corrections to the photoionization rates. For typical 
metagalactic fluxes and IGM densities during reioniza- 
tion, hydrogen and helium are mostly ionized, with frac- 
tions xhi £hii and XHci "C xn e n "C xhcIii- Because 
the first ionization potentials of H I (13.6 eV) and He I 
(24.58 eV) are not greatly disparate, we expect the He III 
and H II ionization fronts to coincide for appropriate val- 
ues of ionizing spectral index. For QSO ionizing spectra 
of power-law form, the He II and H I ionization fronts, 
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Four-panel plot of results from our IGM simulations at z = 2.5 (Smith et al. 2010) showing baryon overdensity (top left), 
ionization fractions of HI and Hell (bottom left and right), and their ratio a'Hdl/^'Hl (top right). Filamentary structures in the cosmic 
web are obvious in both top panels, with a mean value, (^Hell/^Hl) ~ 500, corresponding to r\ = fiHell/ n Hl ~ 40 for primordial He/H = 
y = 0.0823. Much of the Hell in these filaments appears to be photoionized gas at 10 4 ' 3 K, but with some contribution from collisional 
ionization at T > 10 5 K. Original paper with higher resolution available at http://casa.colorado.edu/~mshull/HeII-paper.pdf 



driven by the photon fluxes at 1 ryd and 4 ryd, respec- 
tively, will propagate at the same (photon flux-limited) 
velocity for a critical spectral index, 



ratio is (Fardal et al. 1998; Shull et al. 2004), 



Ocrit 



In 4 



V 1 



^HoIII a H^II ThI 



1.80 



(7) 



corresponding to our adopted He/H abundance, y = 
0.0823. The agreement with the ffijT-observed mean 
spectral index of AGN at 1-2 ryd, (a s ) = 1.76 ± 0.12 
(Telfer et al. 2002), suggests that H II and He III ioniza- 
tion fronts normally propagate together. 

In photoionization equilibrium in regions of high ion- 
ization, ithi *C 1 and cchcII ^ 1) the He II/H I abundance 



N m 



K1.77) 



nmi aity r Heii 



'HI 



Jhi (3 + QHcii) 
^Heii (3 + am) 



rpO.042 

J 4.3 



(8) 



The numerical coefficient has been increased from 1.70 
to 1.77, reflecting an updated value of the primordial 
He/H abundance, y — 0.0823, which we use as an ap- 
proximation to riHein/^Hii- In this equation, , ctj^jj, 
Thi, and Tneii are the case- A recombination rate coeffi- 
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cients and photoionization rates for H I and He II, and 
Jhi and Jneii are the specific intensities of the ioniz- 
ing radiation field at 1 ryd and 4 ryd. The parame- 
ter T4.3 = (T/10 4 - 3 K) is used to scale the temperature 
dependence of the Case-A recombination rates over the 
range T = 10, 000 K to 20,000 K (Osterbrock & Fer- 

land 2006): cJj£ n » (1.36 x 10~ 12 cm 3 s" 1 )^ ' 694 and 

4f « (2.51 x 10- 13 cm 3 s" 1 )^ - 738 . 

Equation (8) is an approximation, which will break 
down when Xfjoii becomes significant. It also does not 
include the effects of collisional ionization of He II, 
which can become important in the shocked filaments of 
shocked gas at T > 10 5 K, the so-called warm-hot inter- 
galactic medium (WHIM). As noted in Figure 6, our sim- 
ulations suggest that some of the He II in shocked IGM 
filaments may be collisionally ionized gas at T > 10 5 K. 
However, the bulk of the He II is photoionized gas at 
T » 10 4 3±0 - 3 K. 

3.4. Troughs of Strong He II Absorption 

As shown in Figures 2 and 4, we have confirmed the 
presence of three long troughs with strong He II absorp- 
tion (thoII > 5) and no detectable flux. These troughs 
have been noted previously (Reimers et al. 1997; Smette 
et al. 2002; Shull et al. 2004), but the very low COS 
backgrounds allow us to better characterize their optical 
depth. We also detect numerous windows of flux trans- 
mission (e.g., at 1138, 1139, 1144, 1154, 1160, 1174, and 
1179 A - see Figure 2). The G130M data provide accu- 
rate boundaries in redshift. From those, we can deter- 
mine their comoving radial sizes, using WMAP-7 cosmo- 
logical parameters: 

1. (Trough 1) z = 2.751 - 2.807 (61 Mpc), 

2. (Trough 2) z = 2.823 - 2.860 (39 Mpc), 

3. (Trough 3) z = 2.868 - 2.892 (25 Mpc). 

These dimensions, all much larger than 10 Mpc, are in- 
compatible with zones of photoionization around typical 
QSOs. They are also unlikely to exist at the epoch fol- 
lowing the overlap of expanding cosmological I-fronts. 
Recent theoretical modeling has studied the statistics 
of He II troughs (McQuinn et al. 2009; McQuinn 2009; 
Furlanetto 2009a,b; Dixon & Furlanetto 2009). These 
studies examine large-scale fluctuations in the He II Lya 
forest transmission, during and after He II reionization. 
A robust result is that long troughs of large optical 
depth, THoii > 4, cannot be explained in models with 
a smoothly evolving ionizing radiation field (Furlanetto 
& Dixon 2010). Furthermore, the observed 25-60 Mpc 
troughs of high He II opacity between 2.751 < z < 2.892 
cannot be accommodated in standard reionization mod- 
els unless He II reionization is delayed until z < 2.7. The 
He II absorption is therefore expected to be patchy at 
z > 2.7, as we have observed. 

It would be useful to conduct surveys (Worseck et al. 
2007) of galaxies and QSOs near the sight line at these 
trough rcdshifts to identify potential sources of ionization 
(or lack thereof). As discussed in Section 3.2, the He II 
opacity greatly exceeds that of H I; the ratio r\ 3> 1 for 
typical conditions in an IGM photoionized by QSOs of 
low space density. Thus, the propagation and overlap of 



I-fronts driven by He II ionizing continua will likely differ 
from those of hydrogen. The greater He II opacity means 
that filaments of the cosmic web still have th ii = 2 — 3 
at z = 2.4 — 2.7 (see Figures 5 and 6). Thus, the topology 
of overlapping He III ionized zones probably differs from 
that of merging H II regions at z = 6. These theoretical 
results are explored in a later paper (Smith et al. 2010) 
using adaptive-mesh cosmological simulations with the 
Enzo code. 

3.5. Fluctuations in He II/H I Ratios. 

Because reionization involves the stochastic propaga- 
tion and overlap of ionization fronts (I-fronts) , one nat- 
urally expects variations in the abundances of H I and 
He II. The metagalactic flux of high-energy (4 ryd con- 
tinuum) radiation varies primarily because of the low 
space density of the dominant ionizing sources (AGN) 
and varying attenuation through the IGM. The varia- 
tions in He II opacity can be especially large during reion- 
ization, because the 4 ryd continuum is strong near AGN 
(inside the I-fronts) but heavily attenuated in the He II 
troughs, far from the AGN and within the gaseous fila- 
ments of the cosmic web. 

Figure 4 shows both He II and H I absorption, over- 
laid in redshift, illustrating the correspondence of He II 
flux-transmission windows with weak spots (low H I ab- 
sorption) in the H I Lya forest. The non-Gaussian wings 
in the COS line-spread function (LSF) is most significant 
for weak absorbers, but we believe this is completely ac- 
counted for in the convolution of the H I data by the 
COS LSF. The resolving power (R = 18,000) of the 
COS LSF is more appropriately the effective Gaussian 
resolving power (Green et al. 2010). The main difference 
in the convolved and unconvolvcd H I data in Figure 4 
is simply a degradation in resolution by a factor of 2.5. 
The COS LSF contributes to a redistribution of flux in 
the narrowest lines, but this is accounted for with the 
non-Gaussian LSF. 

The He II/H I ratio, r\ w 4th ii/thi, tracks the ratio of 
ionizing fluxes at 1 and 4 ryd (see equation [8]). Table 3 
presents the mean and median values of r\ in the G130M 
and G140L data, and over various redshifts. In different 
columns, we present these statistics for the G140L and 
G130M spectral regions, and for the region (2.73 < z < 
2.87) where the two gratings overlap. To calculate rj(z), 
we convolve the high-resolution VLT data with either the 
medium-resolution (G130M) or low-resolution (G140L) 
COS line-spread functions, depending which dataset is 
representing He II at a particular redshift. The COS 
data are then smoothed by seven pixels (~ 1 resolution 
element) and up-sampled to the finer VLT wavelength 
scale. Optical depths th c ii(z) and thi(z) are calculated 
on a pixel-by-pixel basis. 

Figures 7 and 8 show t(z) and rj(z) and demonstrate 
that He II optical depth saturates in certain regions. 
Over the range 2.4 < z < 2.9, 32% of pixels are sat- 
urated, but this fraction depends strongly on redshift 
and spectral resolution. In the low-resolution data at 
z < 2.73, we find that THoii is almost always unsatu- 
rated. At z > 2.7, we see long absorption troughs in 
which 40-45% of the pixels have THoii > 5. We com- 
puted mean and median values of 77 using both saturated 
and unsaturated pixels. As the saturated pixel values for 
He II are entirely in the numerator, our average n values 
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TABLE 3 

Statistics of HeII/HI Column-Density Ratios a 



Statistic 


G140L 
(2.4 < z < 2.73) 
Low- Res 'n 


G130M 
(z > 2.73) 
Medium Res'n 


G140L 
(2.73 < z < 2.87) 
Overlap Region 


G130M 
(2.73 < z < 2.87) 
Overlap Region 


Mean (77) 
Median (77) 
Mean (log 77) 
Median (log 77) 


47 ±42 

qq+56 

M -19 

1.52 ±0.38 
1 r.1+0.43 

0.39 


209 ± 281 

1.90 ±0.68 
1.94 ±0.74 


105 ± 121 
69±^ 

1.83 ±0.41 

1.84 ±0.38 


208 ± 281 

1.89 ±0.68 
1.94 ±0.74 



a We list mean and median values of the ratio 77 = A^eii/A^HI over redshift ranges covered at low resolution (G140L) and medium 
resolution (G130M). The last two columns show r\ statistics for the overlap range, 2.73 < z < 2.87. All statistics are calculated including 
both saturated and unsaturated values of th c ii- Using only unsaturated pixels reduces (77) by 10-20% at 2 > 2.7, with no change at z < 2.7. 




2.4 2.5 2.6 2.7 2.8 2.9 

redshift 



Fig. 7. — Optical depth for Hell (black) and HI (red) and their ratio, 77 ss 47Hell/Tlli versus redshift. Note the change in spectral 
resolution at z < 2.73, the transition between G130M and G140L gratings. The VLT/UVES data have been convolved with the appropriate 
COS line spread function. We only plot values 7jj e ii < 5. Fluctuations at Az 0.01 correspond to proper distances of 10.8 Mpc at z = 2.8. 
For our detectable ranges of optical depth (0.02 < th c ii < 5 and 0.01 < thi < 5.5) we believe r\ can be determined reliably over the range 
5 < 77 < 500 (G130M) and 10 < 77 < 100 (G140L). The mean values and standard deviations are (77) = 47 ± 42 (for G140L range) and 
(77) = 189 ± 276 (for Gl30M range). See full statistics in Table 3. There may be a slight trend of higher 77 toward lower z. 
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2.75 2.80 2.85 2.90 

redshift 



Fig. 8. — Closeup of optical depth for He II (black) and HI (red) and their ratio 7] at redshifts near the QSO, 2.75 < z < 2.90. In the 
higher-resolution G130M data, r\ varies on scales down to Az 0.003. Grey points indicate locations where THeii > 5, so that ij is a 
lower limit. We believe r] can reliably be measured over the range 5 < r) < 500, considering reasonable limits on 0.02 < tjj c h < 5 and 
0.01 < t hi < 5.5. 
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are lower limits. Using only unsaturated pixels reduces 
{rj) systematically by 10-20% at z > 2.7. Correcting for 
larger true values of THoii > 5 would increase the mean 
statistics. There is no change at z < 2.7. 

We confirm the fluctuations seen by FUSE between 
2.73 < z < 2.90 (Kriss et al. 2001; Shull et al. 2004) on 
small scales, Az w 0.003 — 0.01. Figures 7 and 8 show 
that r\ fluctuates on scales down to Az « 0.01 at G140L 
resolution and Az w 0.003 at G130M resolution. For 
reference, the redshift interval Az = 0.01 corresponds 
to a comoving radial distance of 10.8 Mpc at z = 2.8. 
Regions of hot gas (T > 10 K) might appear as high- 
r] systems, in which H I is far more collisionally ionized 
than He II. However, the COS data also exhibit several 
instances of low- 1] absorbers with n « 5 — 20. Since r] ps 
^Heii/THi, the low- 77 systems appear primarily where the 
H I optical depth is large and He II optical depth small. 
If photoionized, these absorbers may be irradiated by 
very hard radiation fields, with large ratios of Jneii/ Jm 
(equation 8), perhaps near an AGN with a particularly 
flat EUV spectral index, a s < 1. 

Several good examples of low-ry absorbers appear at 
z w 2.755, 2.785, 2.830, 2.865, and 2.904. Most of these 
systems have high thi, the exception being the absorber 
at z = 2.865 in which both th c ii and thi are low, and 
their ratio poorly determined. For high-?7 systems, the 
G130M data exhibit regions with n = 200 — 500, typically 
in gas with large THeii but low thi- Extreme excursions 
in i] may arise from noise, and their detectable range 
depends on spectral resolution. We therefore take care 
to assess the range over which we can measure reliable 
optical depths. For ratios within the range r\ = 10-100 
(G140L) and rj = 5-500 (G130M), we believe nearly all 
the fluctuations are real. Unfortunately, strong troughs 
of He II absorption at z > 2.7 A hinder our ability to mea- 
sure THeii for individual systems with the high-resolution 
G130M. Many of the ^-fluctuations appear at z < 2.7, 
observed with the lower-resolution G140L. 

These small-scale ry-variations have been the topic of 
numerous theoretical papers (Zuo 1992; Fardal & Shull 
1993; Fardal et al. 1998; Bolton et al. 2005, 2006). The 
key focus of these papers is on fluctuations in the small 
numbers of AGN within the characteristic attenuation 
length, corresponding to the "r = 1 sphere" at 4 ryd 
photon energy. Because the He II opacity is so much 
larger than that of H I, one expects large variations in 
the He II ionizing flux, JhcII, on comoving scales of 20-40 
Mpc atz« 2.4-2.9 (Fardal et al. 1998; Furlanetto 2009; 
Dixon & Furlanetto 2009). A second important effect is 
that the observed low-redshift AGN exhibit wide vari- 
ations in their EUV spectral indices (Telfer et al. 2002; 
Scott et al. 2004), which produce strong source variations 
in the flux ratio Jn e ii/ Jm- Finally, radiative transfer of 
the metagalactic ionizing radiation reprocesses the ion- 
izing spectra, producing further variations (Fardal et al. 
1998). What is now firmly established by the new COS 
data are the variations in rj on ~ 10 Mpc comoving radial 
distance scales (Az « 0.01). 

3.6. Proximity Effect and Quasar Stromgren Spheres 
It is surprising to note the sharpness of the He II edge 
at 1186.26 A and the lack of any proximity effect near the 
QSO. If ionizing radiation were escaping from the AGN, 
the surrounding IGM should be ionized within 5-10 Mpc 



of the QSO. The He II absorption extends Az « 0.019 
beyond the previous QSO redshift estimate, z om = 2.885 
(Reimers et al. 1997). As discussed in Section 3.1, this 
dilemma is easily resolved by moving the quasar redshift 
to z = 2.904, an offset (1460 km s" 1 at z = 2.9) that cor- 
responds to 19.5 Mpc comoving proper distance. How- 
ever, the two proposed systemic redshifts, zqso = 2.885 
and 2.904, raise several questions. Does the QSO lie 
inside an H II and He III cavity produced by its pho- 
toionizing radiation? Or, is the QSO ionizing radiation 
absorbed by circumnuclear gas? Is there infalling gas 
near the AGN that produces associated absorption with 
Zahs > 2cm- It is worth investigating these possibilities. 

For H I (1 ryd) and He II (4-ryd) ionizing-photon lumi- 
nosities, 5q H ^ and S^ c \ what Stromgren radii will this 
bright QSQ carve out? These distances are obviously 
limiting cases, in ionization equilibrium, whereas one ex- 
pects the ionization zones around the QSO to involve 
time-dependent propagation of ionization fronts. Indeed, 
the lack of a proximity effect could arise from this source 
just starting to burn its way out of the galactic nucleus. 
Given the strength of the z = 2.9 absorber, we will as- 
sume higher gas densities associated with the QSO and 
scale to an overdensity (Sh ~ 100) at z = 2.9. The H and 
He densities are then: n H = (1.13 x 10~ 3 cm~ 3 )(<5 H /100), 
n e w 1.165n/r, and n Hc = (9.30 x 10~ 5 cm- 3 )(5 H /100). 
At T = 10 4 - 3 K, a { nJ u = 9.1 x lO" 13 cm 3 s" 1 and 
af^ = 1.4 x 10~ 13 cm 3 s _1 . A typical (L*) quasar 
produces a total of ~ 10 56 LyC photons s" 1 , but this 
QSO is much brighter (V = 16.1). Extrapolation from 
1200 A down to the He II Lyman limit (889.5 A at 
ZQSO = 2.904) shows that HE 2347-4342 has an ob- 
served specific flux, F\ w 6 x 10~ 15 erg cm~ 2 s _1 A -1 
or F v w 1.6 x 10~ 27 erg cm~ 2 s _1 Hz -1 . From this 
flux and a luminosity distance, dh ~ 7.7 x 10 28 cm at 
z = 2.904, we estimate an ionizing photon luminosity 



of S, 



(Ho) 



vL v j2> ~ 1.6 x 10 57 photons s 1 in the 4- 



ryd continuum. Scaling to S^^ « 10 57 photons s i , we 
estimate a He III Stromgren radius, 



R 



(HoIII) 

Str 



3S*, 
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_47rn e n He a He ii 



! (4 Mpc) 



? (He 



1/3 



"I 1/3 



10 57 



(<V100) - 2 / 3 T 4 °3 255 . 



A similar calculation, using the observed flux, F\ ~ 1.2 x 
10-15 crg cm -2 s -i A" 1 , at the redshifted (3561 A) H I 

edge, gives S^ « vL u /i w 5 x 10 57 photons s _1 and an 
H II Stromgren radius of 



r: 



(HII) 

Str 



35, 



(H) 



47T n e n H Qhi 



: (6 Mpc) 



1/3 



(H) 



5 x 10 57 s^ 1 



1/3 
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-2/3^0.286 
J 4.3 



We note that the ratio of the extrapolated fluxes at the 
respective ionization limits is fairly small. In the QSO 
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rest-frame, F(912 A)/F(228 A) « 5, indicating a fairly 
hard (fiat) power-law spectrum oc v^ 12 . These esti- 
mated Stromgren radii are a substantial fraction of the 
proper distance of 19.5 Mpc between the two proposed 
QSO redshifts, z = 2.904 and z = 2.885. Even though 
the local ionization zone may not extend as far as these 
equilibrium Stromgren radii, the expected proximity ef- 
fect is minimally visible, except for a narrow feature with 
~ 60% flux transmission at z = 2.9055, just below the 
He II edge (Figure 4). We conclude that ionizing radi- 
ation from this QSO is strongly absorbed in the prox- 
imity of the nucleus, possibly by circumnuclear gas or 
by the strong absorber at z w 2.9. The photoioniza- 
tion cross sections for H I and He II above threshold are 
approximately cr H i = (6.30 x 10~ 18 cm 2 )(F/13.6 eV) -3 
and CTHeii = (1-58 x 10~ 18 cm 2 )(F/54.4 cV)~ 3 . Thus, 
to absorb a substantial amount of the ionizing continua 
requires column densities 

N H i = (1.59 x 10 17 cm" 2 ) [E/13.6 eV] 3 t{E) 
NhcII = (6.35 x 10 17 cm" 2 ) [E/bAA cV] 3 t{E) (9) 

for continuum optical depth r at photon energy E. 

For this very bright QSO, the absence of a He II prox- 
imity effect indicates a significant attenuation of the 4 
ryd continuum, with optical depth r > 3 at mean pho- 
ton energies 1.5 times the 54.4 eV threshold, or NhcII > 
6 x 10 18 cm~ 2 . What are the observational consequences 
of such large column densities of ionized gas? It has been 
suggested (Reimers et al. 1997) that O III resonance lines 
at 303.415 A and 305.596 A might be detectable, but we 
find no obvious candidates. It would be useful to search 
for signs of splittings of these lines, which are separated 
by 2.182 A(l + z) w 8.5 A at z w 2.9. Some of the as- 
sociated absorbers have H I anomalously stronger than 
He II, suggesting that the gas is close to the QSO with he- 
lium mostly fully ionized (He III) . More clues may come 
from the strong He II absorbers at z = 2.900 - 2.904 
(1184.7-1186.2 A) which would be closest to the QSO 
if zqso = 2.904. This system has high optical depth in 
both He II and H I (see Figure 4). Very strong metal- 
line absorption is seen in O VI with weaker absorption 
in C IV, C III, and N V (Reimers et al. 1997; Fechner 
et al. 2004). By fitting higher Lyman-series absorbers, 
Fechner et al. (2004) found log N H i = 16.02 ±0.03, while 
the metal-line absorbers gave log Nciv = 13.66 ± 0.03, 
log N NV = 13.42 ± 0.04, and log N OV i > 14.62. They 
proposed that this absorber is exposed to the strongest 
and hardest radiation from the QSO, and therefore may 
shield the other associated absorbers. However, their 
photoionization models gave log NhcII ~ 16.3, which is 
100 times too small to provide the necessary absorption 
(eq. [9]). The H I column is 10 times too small to shield 
the 1 ryd continuum. 

Fox, Bergeron, & Petitjean (2008) obtained newer 
(high S/N) spectra on these metal- line absorbers at 
z = 2.8916, 2.8972, and 2.9041. They found somewhat 
higher column densities of C IV and N V and a slightly 
higher column density, logA^i = 16.29 ± 0.05 in the 
z = 2.9041 absorber. This column is still insufficient to 
shield the QSO's ionizing continuum. We conclude that 
no satisfactory explanation exists for the absence of a 
proximity effect around this very luminous QSO, other 



than the possibility that it has only recently turned on 
(within the last Myr). 

4. CONCLUSIONS AND IMPLICATIONS 

The high throughput and low background of the Cos- 
mic Origins Spectrograph allow us to probe the He II 
reionization epoch toward HE 2347—4342 at redshifts z 
— 2.4-2.9. Because of photoionization and recombina- 
tion rates, He II is a far more abundant species than 
H I, with typical ratios rj = nHeii/'ft-Hi varying from 10- 
200 in the filaments of the Lya forest. With COS, we 
find that reionization of He II to He III is patchy at 
z = 2.7 — 2.9, with 5-10 A troughs of high optical depth 
(THcii > 5) punctuated by narrow windows of flux trans- 
mission. Previous studies of this QSO sight line (Kriss 
et al. 2001; Shull et al. 2004; Zheng et al. 2004) found a 
gradual decrease in He II optical depth at z < 2.7, sug- 
gesting a He II reionization epoch z T = 2.8 ± 0.2. With 
our new COS data, we find a more complicated picture, 
in which He II is slowly reionized at z < 2.7. 
We now summarize the new observational results and 
astrophysical issues: 

1. We have found clear evidence that the epoch of 
He II reionization is delayed below z = 3, as we see 
patchy ionization down to z w 2.7. For this bright 
AGN, high-S/N, low-background observations from 
COS allow us to measure He II absorption, the UV 
continuum of the QSO, and troughs of high-rneii- 
The high spectral resolution of the G130M grating 
probes absorption fluctuations and the He II ion- 
ization edge, and the G140L grating allows us to 
study He II reionization down to z s=s 2.4. 

2. Comparison of COS/G130M and G140L data with 
previous HST and FUSE data confirms the far- 
FUV continuum, absorption troughs, transmission 
windows, and He II optical-depth below the He II 
edge at 1186.26 A. The COS/G140L fluxes agree 
fairly well with FUSE data at z > 2.7 (Shull et al. 
2004), but below 1130 A (z < 2.7) considerable 
differences appear. The COS fluxes are lower than 
those of FUSE inferred by Zheng et al. (2004), and 
our derived He II optical depths are therefore some- 
what higher at 2.4 < z < 2.7. The broad-band 
optical depths (averaged over Az = 0.1) only drop 
below (th c ii) < 1 at z < 2.4. 

3. We observe optical-depth variations in He II and 
H I that imply fluctuations in the He II/H I abun- 
dance ratio on scales Az « 0.003 — 0.01, corre- 
sponding to 4-10 Mpc in comoving distance. Such 
variations likely arise from AGN source variations 
(1-4 ryd spectral indices) and IGM radiative trans- 
fer effects. The number of low-77 systems (mostly 
low THoii and high thi) seen with COS is lower than 
claimed in FUSE observations. 

4. We see intervals of patchy He II flux-transmission, 
as well as three long (5-10 A) intervals of strong 
He II absorption (25-60 Mpc proper distance) be- 
tween z = 2.751 - 2.807, z = 2.823 - 2.860, and 
z = 2.868 — 2.892. These troughs probably arise 
from dense filaments in the baryon distribution and 
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the scarcity of strong photoionizing sources (AGN) 
within 30-50 Mpc. 

5. At z = 2.7-2.9, probes of He II fraction, iEHeii > 
(10~ 2 )(0.1/<!>He)(THeii/5) are now possible at opti- 
cal depths THcii > 5. The He II absorption is 50- 
100 times more sensitive to trace abundances than 
H I at z — 6, where a neutral fraction a;ni = 10~ 4 
saturates the Lya Gunn-Pctcrson absorption. 

What future work is needed to better understand the 
He II reionization? First, we need to assess whether 
the HE 2347-4342 sight line is typical. The COS- 
GTO team will be observing two more He II AGN tar- 
gets, HS1700+6416 and Q0302-003, which may clarify 
whether the patchy reionization at z < 2.8 is common to 
most sight lines. Second, along the current sight line, we 
have identified many regions of partial flux transmission 
and three long troughs of He II absorption, where sur- 
veys of galaxies and AGN will help clarify the sources of 
ionization. Third, infrared spectra of the [O III] A5007, 
4959 narrow emission lines (at 1.955 /J,m and 1.936 /im) 
could confirm the QSO systemic redshift. Pinning down 
zqso is crucial for understanding the nature of the ion- 
ization state of the gas near the QSO, which shows no 
evidence of a proximity effect. On the theoretical front, 
cosmological simulations with time-dependent H and He 
chemistry and I-front radiative transfer should elucidate 
the structure topology of the filaments and voids that we 
observe in the He II Lya-forest and He II troughs. 

Serious studies of He II reionization have been hindered 
by the lack of many AGN targets sufficiently bright to 
obtain far-UV spectra. The spectrographs on FUSE and 
HST have made good progress on three AGN sight lines, 
but newly discovered He II GP targets (Syphers et al. 
2009a,b) are considerably fainter than HE 2347-4342 
or HS1700+6416. New observations of these fainter tar- 
gets will probably measure only broad-band He II optical 
depths, for which the high-resolution COS spectra may 
provide guidance in interpretation. The long-term study 
of the He II reionization process, from z = 2.3 — 3.3, must 
await the construction of a future large-aperture, far-UV 
telescope in space (Shull et al. 1999; Postman et al. 2009). 
It would be especially useful if this mission had high sen- 
sitivity and moderate spectral resolution (R > 20, 000) 
extending from 1130 A down to the 912 A hydrogen Ly- 
man Limit. This would enable the study of patchy He II 
reionization in He II A304 down to redshifts z>2. 
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